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ABSTRACT 

Time-resolved spectroscopy is performed on eight bright, long gamma-ray bursts (GRBs) dominated 
by single emission pulses that were observed with the Fermi Gamma-ray Space Telescope. Fitting the 
prompt radiation of GRBs by empirical spectral forms such as the Band function leads to ambiguous 
conclusions about the physical model for the prompt radiation. Moreover, the Band function is often 
inadequate to fit the data. The GRB spectrum is therefore modeled with two emission components 
consisting of optically thin nonthermal synchrotron radiation from relativistic electrons and, when 
significant, thermal emission from a jet photosphere, which is represented by a blackbody spectrum. 
To produce an acceptable fit, the addition of a blackbody component is required in 5 out of the 8 cases. 
We also find that the low-energy spectral index a is consistent with a synchrotron component with 
a = —0.81 ± 0.1. This value lies between the limiting values of a = —2/3 and a = —3/2 for electrons 
in the slow and fast-cooling regimes, respectively, suggesting ongoing acceleration at the emission 
site. The blackbody component can be more significant when using a physical synchrotron model 
instead of the Band function, illustrating that the Band function does not serve as a good proxy for 
a nonthermal synchrotron emission component. The temperature and characteristic emission-region 
size of the blackbody component are found to, respectively, decrease and increase as power laws with 
time during the prompt phase. In addition, we find that the blackbody and nonthermal components 
have separate temporal behaviors. 

Subject headings: acceleration of particles — gamma-ray bursts — gamma rays: stars — methods: data 
analysis — radiation mechanisms: non-thermal — radiation mechanisms: thermal 
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1. INTRODUCTION 



The fireball model of GRBs (ICavallo fc Rees I [l978l: 



lGoodman1ll986t iPaczvhski 1119861 ) assumes that a large 
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amount of energy is released in a small space leading to a 
fireball that emits 7-rays. The progenitors of these events 
are not known, but are believed to be the collapse of 
massive stars, or the coalescing of two compact objects. 
Details of this model rely heavily on assumptions about 
the initial parameters of the fireball including the radii of 
emission, the baryon load an d the magnetic fi eld associ- 
ated with the plasma outflow (jMeszaros II2006L for review 
see). The observed spe ctra of GRBs are typically inter- 
preted as non-thermal (|Piranlll999t i Zhang et al II2012D 
though there are some exceptions (jRvde 1120041) . This 
interpretation indicates that a dissipative process such 
as internal shocks energizes the electrons to non-thermal 
distributions. Though some versions of the fireball model 
predict that the entire spectrum in the 7-ray band be in 
the form of ther mal emission, we concentr ate on the ver- 
sion presented in iMeszaros fc Reesi ()2000[ ) where a mix- 
ture of thermal and non-thermal (synchrotron) emission 
emerges. Alternative models such as the Poynting-flux 
dominated outflow (PFD) model exist and could provide 
a viable mechanism for generating the observed prompt 
emission. These PFD models have not yet advanced to 
quantitative spectral and temporal predictions; there- 
fore, it is difficult to compare our results to theses models 
as extensively as is possible for the internal shock model 
for which extensive simul ations of lightcurves and spec- 
tra have been developed (iDaigne fc Mochkovitch 1 119981 : 
iBosniak. Daigne. fc Dubus II2009D . 

The physical emission process behind these extreme 
cosmic explosions has not been established. Histori- 
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cally, GRB spectra have been best fit by a parametrized smoothly broken power-law that is connected expon en- 

tially, known as the Band function (iBand et al.lll993[) 
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where £ is the photon energy. The Band function is 
characterized by its low and high-energy power-law in- 
dices (Band a and /3 respectively) as well as its vF^, 
peak energy, Epcak- The ability of the Band function 
to fit most prompt-spectra has led to the ubiquitous use 
of the Band function's fitted parameters as indicators of 
the underlying emission and acceleration processes. In 
particular, the value of the low-energy index, Band a, 
has been of interest because non-thermal spectra differ 
most noticeably below their respective vF^ peaks where 
the spectral index is fixed based on the type of emis- 
sion (synchrotron, inverse-Compton, etc). Above the 
peak, the high-energy index is variable and related to 
the distribution of the emitting electrons and provides 
no clues for discerning between emission processes. Al- 
thoug h such extrapolations can be useful, iBurgess et al. I 
(|2011[ ) (hereafter 1311) directly fit a synchrotron photon 
model to the prompt emission data of GRB 090820A, 
showing that it is possible to fit a physical model di- 
rectly to the data without relying on interpretations of 
the Band function in order to determine the process of 

7-ray emission. 

P re yious studies (iGonzalez et al. ll2003HGuiriec et ahl 



IMnl : IGuiriec et al.ll2nTl IRvde et alJl2ninHRvde 1120051) 
showed that a single Band function cannot fully account 
for the spectrum of all GRBs. The addition of a black- 
body and/or power-law better describes the spectra of 
these GRBs. Specifically, the addition of the blackbody 
below Epcak also allows for the direct fitting of the syn- 
chrotron model to the data, which we exploit in this 
work. The addition of a blackbody in the spectrum al- 
lows a quantitative analysis of its properties to ascer- 
tain its origin. Several studie s JAIeszaros & Rees 200(1 
Daigne fc Mochkovitch I [20021 iRvde 1 120041 IRvde et alj 



2008f ) have developed the theoretical framework for a 



blackbody component coexisting with a non-thermal 
component in the spectra of GRBs, as well as testable 
relations that can be applied to data to investigate the 
photosphere of GRBs. 

Herein, we extend the analysis of Bll to several bright 
Fermi GRBs and in addition we investigate the ability of 
the Band function to serve as a proxy for physical emis- 
sivities in the fitting process. We will show that slow- 
cooled synchrotron is a viable model for GRB prompt 
emission using data from both the Gamma-ray Burst 
Monitor (GBM) (10 keV to 40 MeV) and Large Area 
Telescope (LAT) (100 MeV - 300 GeV) on board the 
Fermi space telescope. This analysis makes use of the 
newly released LAT low-energy (LLE) data which pro- 
vides a larger effective between 30 MeV and ~1 GeV al- 
lowing us to extend the LAT analysis down to 30 MeV. 
First, we give a description of the non-thermal and black- 
body photon emissivities (Section [2|) . We then present 



our analysis technique and compare the results with pre- 
dictions about each component (Sections [3] and 0]). 

2. MODEL SPECTRAL COMPONENTS 

In the fireball model of GRBs, the majority of the 
flux is theoretically expected to be in the form of 
thermal emission coming from the photosphere of the 
jet. However, over 2/3 of the low-energy indices re- 
cover ed by the Band functio n arc too soft to be ther- 
mal (jGoldstein et al. ]|2012ai rbl). This points to a non- 
thermal emission process for most GRBs. Multi-spectral 
component analysis of Fermi GRBs has shown that 
while the majority of the emission is non-thermal, a 
small fraction of the energy radiat ed apparently origi- 
nates from a blackbod y component (jGuiriec et al.ll201l] : 
lAxelsson et al. Il2012|) . In Bll, a blackbody component 
was also identified when the non-thermal emission was 
fit with a synchrotron photon model. This combination 
of blackbody and non-th ermal emission was predicted by 
iMeszaros fc Reesi ()2000D . We further improve the anal- 
ysis of Bll by first reviewing the synchrotron emission 
from cooled electrons and the slow-cooled synchrotron 
model from Bll. In addition, we review several ob- 
servable relations of the blackbody component. These 
components are then implemented into a fitting program 
which directly convolves the physical models with the 
GBM detector response to compare against the observa- 
tions. 

2.L Synchrotron Radiation 

While no definitive model for the non-thermal emis- 
sion of GRBs exists, synchrotron is the simplest and 
most efficient non-thermal emission process that could 
account for the observed GRB radiation. It is a natu- 
ral choice in the evolution from empirical fitting func- 
tions to physical models. In Bll, we implemented a 
parametrized synchrotron model for fitting GRB data 
directly. We implemented a synchrotron model first 
used in fittin g the deconvolve d spectra of GRBs b y 
Tavanil (I1996D and then later b y iBaring fc Brabvl (|2004D . 
Sari. Piran. fc Naravan I ()1998[ ) identified two regimes for 
the synchrotron emission which arc relevant for GRBs: 
fast and slow-cooling. The difference between the two is 
related to the radiative time-scale of the emission. 

For the slow-cooling synchrotron mod el, we assume the 
electro n distribution from Bll (see also IBaring fc Brabvl 
(|2004[ )) which includes electrons in a thermal pool as well 
as electrons that are accelerated into a power-law tail (see 
Figure [l|) . 
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Here, no is the normalization, 7 is the electron Lorentz 
factor, 7th is the thermal electron Lorentz factor, 7„iin 
is the minimum electron Lorentz factor of the power- 
law tail, e is the normalization of the power-law, S is 
the electron spectral index. The function 0{x) is a step 
function where Q{x) = for a; < 1 and 8(x) = 1 for 
a; > 1. There is strong evidence that sub-relativistic 
shock acce leration produces similar electron distributions 
to eq. [2] (|Baring. Ellison. &: Jones 1 11995() . Recent de- 
velopments in relativistic shock acceleration have shown 
that it is possible to accelerate elect rons from a ther- 
mal pool into a power-law distribution (ISpitkovskv1l2008l : 
iBaring |[20lTt iSummerlin fc Baring Il2012[ ). Theseresults 
motivate our choice of electron distribution even though 
our parametrization may oversimplify the process of ac- 
celeration and ignores the effects of photon-electron col- 
lisions such as Compton scattering that are likely taking 
place in the high 7-ray regime. Yet, the distribution we 
chose closely resembles the expected distribution result- 
ing from shock acceleration in the absence of synchrotron 
cooling. The inclusion of the thermal electron component 
is physically required and is much closer to reality than 
the simple power-law electron distribution commonly in- 
voked in GRB modeling. Simulations that include the 
full radiative heating and cooling of electrons expected in 
the high 7-ray regime but that neglect acceleration pro- 
cesses produce electron distributions that differ greatly 
from simple power-laws (jPe'er fc Waxmaini2004D . These 
complicated electron distributions are difficult to model 
in a fitting process and may require spectral resolution 
beyond what the Fermi data provide to identify in spec- 
tral fits. 

We convolve this simplified distri bution with the stan- 
dard isotropic synchrotron kernel (jRvbicki fc Lightmanl 
[1971 
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expresses the single-particle synchrotron emissivity (i.e., 
energy per unit time per unit volume) in dimensionless 
functional form. The Kc^i-^ term is the modified spheri- 
cal Bessel function. The characteristic scale for the syn- 
chrotron photon energy is 
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where E'* = |-g^TOeC^, B is the global magnetic field 
strength, F is the bulk Lorentz factor, and Bex = 4.41 x 
10^'^ Gauss is the quantum critical field. 

In principle, there are six spectral parameters that can 
be constrained by the fits: no, i?*, 5, e, 7th, and 7niin; 
however, we fix 7th, 7min, and e due to fitting correlations 
as explained below. The parameter E^ scales the energy 
of the fit and is linearly related to the Band function's 
Epeak- Numerical simulations of GRB s hocks have shown 
that the ratio of 7th and 7min is ~3 ([Baring &: Brabvl 
|2004[) and therefore we fix those parameters to have that 
ratio. The parameters i?,, 7th, and 7,nin all directly scale 
the peak energy of the spectrum but do not alter its shape 
and thus cannot be independently determined. For this 



reason we chose values of 7th = 10 and 7niin — 30 for 
all fits and left Et, free. Estimates of B, the magnetic 
field, and F are also not obtainable from these spectral 
fits because of the correlation between _E* and the elec- 
tron Lorentz factors. For this reason the values of £'* 
from our fits are useful only for examining the evolution 
of the peak energy and not for obtaining the magnitude 
of B and F. It was found in Bll that the parameter 
e, corresponding to the relative amplitude between the 
thermal and non-thermal portions of the electron dis- 
tribution, was not easily constrained in the fitting pro- 
cess and produ ced non-physical elec tron distributions as 
pointed out in iBeloborodov I (|2012[ ). Therefore, we nu- 
merically fix this parameter so that there is no disconti- 
nuity between the thermal and non-thermal parts of the 
distribution. After these simplifications, three shape pa- 
rameters remain free: E', , 5 and no, which corresponds to 
the amplitude. Compared with the Band function's four 
fit parameters this model is simpler yet tied to actual 
physical processes. 

The second regime of synchrotron emission we exam- 
ine is the so-called fast-cooling regime. Synchrotron is 
a very efficient process and cools electrons rapidly. In 
fact, the dynamic timescales associated with GRBs are 
much lon ger than the cooling timescale s of synchrotron 
emission ([Sari. Piran. fc Naravanlll998[) . It is therefore 
necessary to investigate synchrotron emission from elec- 
trons that have been rapidly cooled. We assume that 
some acceleration process injects a power-law distribu- 
tion of electrons 
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into a region where they are allowed to cool. Here, qe is 
the electron density. We neglect the inclusion of the ther- 
mal pool because its association with the radiative region 
is poorly understood. Moreover, we will show that the 
fast-cooling spectrum is already too broad for the typical 
GRB vFi, peak and the thermal distribution only broad- 
ens the spectrum further. The coo fing of electrons is gov- 
erned by the continuity equation ([Blumenthal fc Gould I 
[T970I) 
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where Ue/tesc represents the loss of particles from the 
emission region from which we can define the maximal 
cooling scale 7/7 ~ tesc- This corresponds to a Lorentz 
factor, 7coo;, below which cooling shuts off. We expect 
that 7coo; lies well below ^rain owing to the fact that a 
cooling break has not been observe d in GRB data (how- 
ever, see lAckermann et al. I (|2012bD ). For very high elec- 
tron Lorentz factors, 7/7 ^ t^sc and therefore the loss 
term can be neglected, i.e., the dynamical timescale is 
much longer than the radiative timescale. With this as- 
sumption we can simplify eq. [7| and it becomes time in- 
dependent. The resulting electron distribution can then 
easily be solved for 
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Substituting in the synchrotron coohng rate 
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where r^ — rrieC^ / {eB) and tq — e^ jim^c')^ eq. [5] 
yields the synchrotron-cooled broken power-law distri- 
bution of electrons (see Figure [2|) 
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This distribution is convolved with eq. H] to produce a 
photon spectrum. The radiation spectrum generated by 
this distribution has an asymptotic low-energy index of 
—3/2 , the so-called 'second line-of-death' (i Cohen et aT~1 
I1997D . The free parameters used for fitting this spectrum 
are the electron spectral index, (5, the fast-cooling equiv- 
alent of Et, and the overall amplitude. However, the ma- 
jority of GRB spe ctra in the BATSE (iGoldstein et al. I 
I2012b[) and GBM (jGoldstein et al. Il2012al) spectral cat- 
alogs have spectra with low-energy indices harder than 
—3/2. Even if the spectra were consistent with the 
low-energy index, the fast-cooling spectrum's curvature 
around the vFy peak is much broader than the peak of 
the data (See|lJ|). 

Such simplified models cannot capture all the nuances 
of a detailed simulation of the acceleration and radiation 
processes in GRBs. Both the cooled and uncooled 
models neglect self-absorption and Klein-Nishina effects 
which may be significant but which are beyond the scope 
of this investigation. One possible alternative to fitting 
parameterizied models is forming templates of spectra 
generated by a simulation and fitting those to data. 
The models implemented here already have degeneracies 
that require the fixing of certain parameters due to 
the lack of spectral resolution in the data. Simulations 
typically have more param et ers than our current model 
dPe'er fc Waxman ' 2004; 'Asano. Ino ue. fc Meszarosi 
[2009; Bosniak, Daigne. & Dubus 2009) and constrain- 
ing those models via spectral templates using data from 
Fermi will be difficult. 

2.2. Blackbody Component 

The pure fireball scenario for GRB emission pre- 
dicts that most of the flux is fro m thermal emission 
(| Goodman I II 9861 ; IPaczvhskTI 119860 . This is because as 
the jet becomes optically thin at some photospheric ra- 
dius, rph, it releases radiation that has undergone many 
scatterings with the optically thick electrons below the 
photosphere. We model this emission as a blackbody 
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where A is the normalization and kT scales the energy 
of the function. This is simplified thermal emission 
from the photosphere that does not take into account 
the effects of relativistic broadening that can produce 
a multi-color blackbody emerging from t he photosphere 
(IBeloborodov I [2010l: IRvde et al.l IMot IPe'er fc Rvdel 
120 111) . iRvde et al. I (I2008D showed that if it is assumed 
that the thermal component is emanating from the pho- 
tospheric radius of the jet, several properties about the 



blackbody component are derivable. In particular, the 
ratio of the thermal [Ubb) and non-thermal {Unon-th) 
flux is given by 
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where p is the density of the plasma, cr^f, is the Stefan- 
Boltzmann constant, and kT oc p^/'^. The cooling be- 
havior is also well predicted for a thermal component. 
The temperature should decay as T ex t^"^/^ . It has been 
found that the evolution of kT often follows a broken 
power-law trend with the index below the break averag- 
ing to 2/3 (|Rvde and Pe'er II2009D . Finally, a true 

blackbody has a well defined relation between energy flux 
and temperature: 



Fbb = Nasbr 



(13) 



where iV is a normalization accounting for geometrical 
effects. 

The photospheric radius and the transverse size of the 
photospheric emitting region are also of great importance 
to understanding the geom etry and energetics of GRBs. 
In lRvde and Pe'erl ()2009[ ). a parameter TZ 
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is used to track the outflow dynamics of the burst. Sev- 
eral BATSE GRBs were found to have a power-law in- 
crease of TZ with time. However, the connection between 
TZ, T, and Fbb was difficult to establish because the er- 
ror on the data points was large. Understanding these 
connections is essential to unmasking the structure and 
temporal evolution of GRB jets. 

3. TIME RESOLVED ANALYSIS 

3.1. Summary of Technique 

The GRBs in our sample were selected based on 
two criteria: large peak fiux and single-peaked, non- 
overlapping temporal structure. The GRBs were binned 
temporally in an objective way described in Section 13.31 
and spectral fits were performed on each time bin using 
four different photon models (Band, Band-fblackbody, 
synchrotron, synchrotron-|- blackbody). When fitting 
synchrotron we compared the fits of slow-cooling and 
fast-cooling synchrotron. In many cases the fits from 
fast-cooling synchrotron completely failed. From these 
spectral fits a photon fiux lightcurve was generated for 
each component and fitted with a pulse model to deter- 
mine the decay phase of the pulse. We describe each step 
in the following subsections. 

3.2. Sample Selection 

To fully constrain the parameters of the fitted models, 
we selected GRBs with a requirement that the peak flux 
be greater than 5 photons s^^ cm~^ between 10 keV and 
40 MeV. For the GBM analysis, the time-tagged event 
(TTE) data was used due to its high temporal and spec- 
tral resolution. It is important for our GRBs to have 
a simple, single-peaked lightcurve structure to avoid the 
overlapping of different emission episodes. While we can- 
not be sure that a weaker emission episode does not lie 
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beneath the main peak, the bursts we selected have no 
significant additional peak during the rise or decay phase 
of the pulse. These two cuts left us a sample of eight 
GRBs: G RB 081110A (iMcBreen et ah I [2008I ). GRB 
081224A (j Wilson-HodgeejaLJ [2QQ|)7~GRB 090719A 
(Ivan der Ho rst 2009). GRB 090809B, GRB 100707A 
(jWilson-Ho dge et all [20Tol). GRB 110407A, GRB 
110721A, f,Ticrncv et al. I |201H) and, GRB 110920A 
(Figure [g. GRB 081224A and GRB 110721A were 
both analyzed including the new LAT Low-Energy (LLE) 
data that provides a high effective area above 30 MeV 
for the analysis of short-lived phenomena, thanks to 
a loose ned set of cut s with respect to LAT st andard 
classes (jPelassa 1 120091 lAckermann et al. |[2"012aD . This 
data selection bypasses th e typical photon classification 
(jAckermann et aLll2012cl ) tree and includes events that 
would normally be excluded but can be selected tempo- 
rally when the signal to background rate is high, such as 
with GRBs. GRB 081224A had very little data above 
30 MeV but the LLE data helped to constrain the spec- 
tral fits. From this sample, five GRBs (GRB 081224A, 
GRB 090719A, GRB 100707A, GRB 110721A, and GRB 
110920A) had blackbody components that were bright 
enough to analyze (Tabled]). 

3.3. Time Binning 

In order to bin the count data, the time bins for spec- 
tral an alysis of each GRB were chosen using Bayesian 
blocks (jScargle et al.l 120131 ). Bayesian blocks define the 
time bins by looking for significant changes in the data 
count rate that define change points based on a Bayesian 
prior. For GBM data, we combined the TTE data from 
the brightest Nal detector with the TTE data from the 
BGO detector and ran a Bayesian block algorithm to find 
the times of the change points. We found that a prior of 
8 gave a good balance between time resolution and hav- 
ing enough counts to perform spectral analysis. These 
change points were mapped to the other detectors used 
in the analysis. This method only accounts for changes 
in counts and therefore could inadvertently combine bins 
with different spectral shapes. 

3.4. Spectral Analysis 

For spectral fits we used the RMFIT ver4.]Ej software 
package developed by the GBM team. Fitting the syn- 
chrotron photon model requires a custom module de- 
veloped and used in Bll. Each time bin was fit with 
one of the four spectral models mentioned above. An 
off pulse region was selected as background and fit with 
a polynomial in time for subtraction for both the TTE 
and LLE data. We fit the physical models to compare 
the validity of each one against the other and the Band 
function to try and understand how the Band function 
parameters correlate with the best fit physical model. 
If the addition of blackbody component did not make 
a significant improvement of at least 10 units of C-stat 
(jArnaud. Smith, fc Siemiginowskal 120 lit ) for any time 
bins of a particular GRB, then we did not include the 
blackbody component in the analyzed fits for that burst. 
C-stat is a likelihood statistic that is modified by a con- 
stant that causes the statistic divided by the degrees of 

^"^ http://fcrmi.gsfc.nasa.gov/ssc/data/analysis/uscr/ 



freedom to approach a x^ value to the limit of a large 
number of counts. However, near the end of the prompt 
emission in some GRBs, the blackbody component be- 
comes weak but has spectral evolution consistent with 
more significant time bins in the burst. The spectral pa- 
rameters of the blackbody in those bins were included 
even though they contributed large error bars to some 
quantities. We checked with simulations that this cut 
was sufficient to identify a significant addition of a black- 
body to the fit model. 

4. RESULTS 

4.1. Test of Slow- Cooling Synchrotron 

In nearly all cases the synchrotron or syn- 
chrotron+blackbody model produced a fit with a 
comparable or better C-stat than the Band function. 
The GRBs exhibited hard to soft spectral evolution 
(Figure H]) for both components. From these fits we 
can derive several interesting properties of the bursts. 
The results of fitting the non-thermal part of each 
time bin in our sample with slow-cooled synchrotron 
indicate that this model can indeed fit the data well. 
The C-stat fit statistic per degree of freedom was at or 
near 1 for most time bins. The spectral fit residuals 
cluster around zero, with no deviations at low-energy 
that might indicate an additional power-law component 
(Figure [5]). The residuals are below Aa for the entire 
energy range. As an example, the fit C-stats for GRB 
100707A and GRB 110721A are shown in Tables H 
and [3] respectively. They show that the slow-cooled 
synchrotron model fits the data as well as the Band 
function when a blackbody is included. The fit C-stats 
for fast-cooled synchrotron are shown for GRB 110721A 
to compare all three non-thermal models. These results 
imply that slow-cooled synchrotron is a viable model for 
GRB prompt emission. We cannot claim it provides a 
better fit to the data than other untested models and 
we will investigate and compare other physical models 
in future work 

An important parameter constrained in these fits is 
the electron index, 5, of the accelerated power-law. The 
canoni cal value from relativistic shock accelerat ion is 
6=2.2 (|Kirk fc Schneider |[l987l: iKirk et al. l[2000l) . The 
distribution of constrained 5's (Figure [6]) is broad and 
centered around (5=5. This steep index could provide 
clues for the structu re and magnet ic turbulence spec- 
trum of the shocks. iBaring I (|2006| ) shows that shock 
speed, obliquity, and turbulence all have a strong effect 
on the electron spectral index of the accelerated elec- 
trons. Steeper indices correspond to increasing shock 
obliquity. Fit models which are built from the electron 
distribution such as the one used in this work enable a 
direct diagnostic of the GRB shock structure. 

4.2. Test of Fast- Cooling Synchrotron 

In order to see if any spectra were consistent with 
the fast-cooling synchrotron spectrum, we implemented 
a fast-cooled synchrotron model where the electrons were 
distributed according to the broken power-law in eq. 1101 
As with the slow-cooling fits, 7min was held fixed to 30. 
Several spectra were tested and all resulted in very poor 
fits regardless of whether the low-energy index found 
with the Band function was much harder than —3/2 (Fig- 
ure [7] and Tabic [4]). This is due to the broad spectral 
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curvature of the fast-cooled spectrum around the vFi, 
peak. The broken power-law nature of the electron dis- 
tribution is smeared out by the synchrotron kernel and 
cannot fit the typical curvature of the GBM data. In 
fact, the fast-cooling synchrotron spectrum has a spec- 
tral index of —2/3 below the 7c which we have fixed at 
1. The fitting algorithm increased E^ to high values to 
align the —2/3 index with the data, which resulted in 
poor fits (Figure HI) . Even when a blackbody is present 
in the bursts, fast-cooled synchrotron is not a good fit 
to the non-thermal part of the spectrum (Table [3]). The 
lack of GRBs with low-energy indices as steep as —3/2, 
additionally disfavors fast-cooled synchrotron as the non- 
thermal emission component in GRB spectra. 

4.3. Synchrotron vs. Band 

The Band function has been used in the literature 
as a proxy for distinguishing among non-thermal emis- 
sion mechanisms. The predicted non-thermal emission 
of GRBs is typically characterized as a smoothly-broken 
power-law with the high-energy spectral index related 
to the index of accelerated electrons and the low-energy 
index related to the radiative emission process. There- 
fore, fitting a Band function to the emission spectrum 
of a GRB should serve as a diagno stic of the radiative 
process responsible for the emission. iPreece et al.l (|1998[ ) 
examined the BATSE GRB catalog and looked at the 
distribution low-energy indices from Band function fits. 
They found that the distribution peaked at a « — 1 and 
that 1/3 of the fitted spectra had low-energy indices too 
hard for synchrotron radiation. The assumption is that 
the Band function's shape approximates synchrotron but 
has an added degree of freedom in the low-energy in- 
dex. However, the Band function has a broader range 
of curvatures around the i^F,^ peak allowing it the pos- 
sibility to deviate from the shape of synchrotron above 
and around the vF^, peak. The synchrotron vF^, peak 
is ex 7min^ 0^ ^*> leading to the relation between Band 
and synchrotron models Epcak cc E^ . This relationship is 
easily recovered from our sample (Figure |9|). Direct com- 
parison of the quality of the fits using Band and the syn- 
chrotron model is not the goal of this study. Both Band 
and the synchrotron model fit the data well with their 
respective fit residuals not deviating more than 4tT and 
centered around zero (Figure [5]). It is important to stress 
that the questions being asked arc does the synchrotron 
model fit the data? and, what temporal evolution do the 
synchrotron parameters undergo? 

For all GRBs in our sample that include both a black- 
body and non-thermal component we compare the pho- 
ton flux (photons s~^ cm~^) lightcurves (integrated from 
10 keV - 40 MeV) derived from synchrotron fits with 
those derived from Band fits (Figure [TU)) . It is seen 
that while both methods recover the same total flux, the 
flux from the individual components is much better con- 
strained when using the synchrotron model for the non- 
thermal component. This is due to the pliability of the 
Band function below Epcak that is not afforded to the 
synchrotron model. 

The C-stat fit values for the synchrotron model loosely 
correlate with the value of Band a found by fitting the 
same interval with the Band function. When Band a was 
much harder than zero, the synchrotron fit was poor and 
typically required adding blackbody to fit the data. The 



flexibility of the Band function with its low-energy power 
law creates the possibility that the index alpha of that 
power law will not accurately measure the true slope if 
Epoak is too close to the low-energy boundary of GBM 
data. Simulated spectra using the Band function were 
created with a grid in both Band a and Epeak to ensure 
that low values of Epeak do not affect the reconstruction 
of Band a in our fits. It was found that Band a could 
be accurately measured when Epeak was as low as ~20 
keV. While the asymptotic value of synchrotron is —2/3, 
fitting the photon model with an empirical function like 
Band with a slightly different curvature could result in 
measured low-energy indices that are different. To mea- 
sure this effect, simulated synchrotron spectra with dif- 
ferent E^, were fit with the Band function. The Band 
a showed a slight dependence on the synchrotron peak; 
moving to softer values for lower E^. The distribution 
of fitted Band a values from these simulations centered 
around —0.81 ± 0.1, a slightly softer value than —2/3 
which may explain the clust ering of Band a at —0 .82 
in the GBM spectral catalog (jGoldstein et al. Il2012"a[) if 
a majority of the non-thermal spectra are the result of 
synchrotron emission. 

4.4. High-Energy Correlations 

There is a well-known spectral evolution in GRB 
pulses of Epcak evolving from hard to soft (see Fig- 
ures [3] and 13]). This leads to two time- resolved cor- 
relations between hardness ( measured as Eneak ) and 
flux (iGolenetski et all 119831: iLiang fc KargatisT 119961: 



iGhisellini et al. 1120101 ). fliang fc Kargatis I (|1996r ( here- 
after LK96) showed that the hardness intensity correla- 
tion (HIC) which relates the instantaneous energy flux 
Fe to spectral hardness and can be deflned as 



Fe — Fn 



E, 



peak 



E, 



'peak,0 



(15) 



where Fq and Epcak, o s-re the initial values at the start 
of the pulse decay phas e and q is the HIC index. 
IBorgonovo fc Rvde I (|2001D found that 57% of a sample 
of 82 BATSE GRBs were consistent with this relation. 
The second relation is the hardness-fiuence correlation 
(HFC) which relates hardness to the time-running flu- 
ence of the GRB. Time-running fiuencc, $(t), is defined 
as the cumulative, time-integrated flux of each time bin 
in a GRB. The HFC is expressed as 



-^pcak -^pcak.O^ 



-*(t)/*o 



(16) 



where $o is the decay constant. LK96 noted that this 
equation is similar to the form of a confined radiating 
plasma. This should not be the case for optically-thin 
synchrotron. Upon differentiating cq. [TH]it becomes ap- 
parent that the change in hardness is nearly equal to the 
energy density: 
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The HFC could be the result of a confined plasma with 
a fixed number of particle cooling via 7-radiation as pro- 
posed by LK96. Since these relations are only valid dur- 
ing the decay phase of a pulse the value of T^ax (the 
time of the peak flux) from the pulse fit of each GRB is 
used as the initial point for Fq and i?peak,o- 
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The use of a hardness indicator is somewhat ambigu- 
ous. Historically, the ratio of counts in low and high- 
energy channels was used as a hardness measure. This 
has an advantage of being model-independent but suf- 
fers from the lack of information associated with the in- 
strument response. High-energy photons can scatter in 
the detector and not deposit their full energy thereby 
artificially lowering the hardness ratio. LK96 used the 
Band function Epeak to compute hardness, which as a 
deconvolved quantity is less instrument-dependent but 
introduces a model dependence. We take this approach 
for both the Band and synchrotron model fluxes. For 
synchrotron we use the _E* parameter as our hardness 
indicator. This is justified by the relationship between 
E^ and Epoak (see Sec. 14.31 and Figure [9|). 

We compute the HIC and HFC for the synchrotron fits 
for each GRB in our sample (Figure [Til and Table [5]). All 
the GRBs seemed to follow the HIC to some extent. We 
find that the HIC index for E^, ranges between «l-2, with 
1 being expected if the emission is due to synchrotron. 
When using the Band function Epeak as a hardness in- 
dicator, it is expected tha t Epeak oc L^/^ oc Fp' from 
basic synchrotron theory (jGhisellini et al. ||2010D . The 
synchrotron fits seem to obey the HFC fairly well. Owing 
to the large errors in the Band flux, fits with synchrotron 
are more consistent with the HFC and HIC than those 
with Band. The deviations of the data from the expected 
synchrotron HIC may be due to the fact that there are 
overlapping pulses under the main emission that alter the 
decay profile. In addition, the use of Bayesian blocks to 
select time bins ignores spectral evolution. If bins with 
very different Epeak are combined then it could affect the 
the HIC and HFC data. 

4.5. Blackbody component 

For most of the spectra in our sample, the black- 
body's iyF,y peak is below the i/F^ peak of the non- 
thermal component. There is sometimes a much 
larger change in C-stat between fits with synchrotron 
and synchrotron-|-blackbody than those of Band and 
Band+blackbody owing to the fact that the Band func- 
tion has more freedom in the shape below Epeak (Ta- 
bles ini and [7]). Simulations of both Band and slow- 
cooling synchrotron were used to find the significance 
of adding a blackbody to the spectrum. It was found 
that even when the difference in C-stat between Band 
and Band+blackbody was greater than the difference be- 
tween synchrotron and synchrotron-f-blackbody, the sta- 
tistical significance in the goodness of fit after the addi- 
tion of the blackbody as high if not greater for the syn- 
chrotron-t-blackbody model for many cases. Computa- 
tional time limits kept us from checking if the significance 
reached 5a. We now focus on the blackbody component 
that is found in the synchrotron-|-blackbody fits. The 
blackbody appears to have a separate temporal struc- 
ture from the non-thermal component, typically peaking 
earlier in time and decaying before the non-thermal emis- 
sion. (Figure [T0|). 

The form of the blackbody used in this work (eq. 
Ilip is simplified and therefore will likely only approx- 
imate the true form of thermal emission from a GRB 
photosphere. Since the blackbody is weaker than the 
non-thermal (synchrotron) component in the spectrum. 



the effects of a broadened and more realistic relativis- 
tic blackbody are masked and only slightly affect the fit 
when combined with the synchrotron model. However, 
in the c ase of G RB 100707A, the blackbody is very bright 
(Figure 10(c) ) and subtle changes in actual shape of the 
photospheric emission become more apparent. This is 
refiected in the C-stat values in Table [21 The Band func- 
tion combined with the standard blackbody is a better 
fit than when using synchrotron as the non-thermal com- 
ponent. In this case, the Band function a is still very 
hard indicating that the Band function is making up for 
additional flux that the blackbody is not taking into ac- 
count. To test this hypothesis, we flt the synchrotron 
model with an exponentially cutoff power-law to mimic 
a modified blackbody. We found that the fits were as 
good as the Band function combined with blackbody fits 
indicating that when the blackbody is bright compared 
to the non-thermal emission a more detailed model of the 
photospheric emission is needed to fit the thermal part 
of the spectrum. 

For each time interval the energy flux (ergs s^^ cm^^) 
ratio is calculated for the synchrotron and blackbody 
components and the sum of these is calculated to ob- 
tain the total flux ratio (see Table [8]) . We flt the ratio 
(^bb/^synch)(t) vs kT(t) to sce if it foUows the trend of 
eq. [121 None of the components in our subset followed 
a linear trend but the values of the slope were all uncon- 
strained due to the high error in the energy fluxes. The 
HIC for the blackbody component is equivalent to eq. [T3l 
however, nea rly all the blackbodies had a HIC index o f 
q K, 2 (Figure [l2(a)|). These results conflrm|^de] (|2004[ ) 
and lRvde and~Pe^r I (|2009l ) who fit BATSE spectra with 
a combination of a blackbody an d a power-law to acc ount 
for the non-thermal component. iRvde et al. I (|2008[ ) de- 
scribes a toy model for the blackbody that allows a range 
of temperature indices related to the internal structure 
of the photosphere that may account for these results. 

Another interesting quantity that can be obtained from 
the blackbody is the HFC. All GRBs in our blackbody 
subset had blackbodies consistent with the HFC (Figure 
112(b)). The decay constants were all of similar value. 
ICrider et al. I (|1999( ) noted that similar values of $o for 
non-thermal components arise as a consequence of a nar- 
row parent distribution. A deeper investigation of a 
larger sample is required to assess if the same is true 
for the blackbody components. 

The temporal evolution of kT for the blackbody of each 
burst appears to follow a broken power-law fFigu re | 13l). 
The e volution is fit with the function derived in iRvde I 
(|2004r ) where we fixed the curvature parameter, (5, to 
0.15. The coarse time binning derived from Bayesian 
blocks does not allow for the decay indices to be con- 
strained for all the bursts but a small subset are close to 
—2/3, as expected (see Table [9]). The temporal decay of 
the blackbody is different than the power-law decay of 
E^^,, indicating a different emission component. 

The TZ parameter was observed to increase with time 
for all the GRBs ( Figure [T4|) . There are breaks and 
plateaus in the trends that do not seem to correlate with 
the breaks observed in the evol ution of kT or w i th the 
flux history of the blackbody. IRvde and Pe'er I ()2OO90 
found that the evolution of TZ can be quite complex but 
mostly follows an increasing power-law that seems inde- 
pendent of the flux history even for very complex GRBs. 
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For those complex GRB lightcurves, it was found that 
analyzing different intervals of the overlapping pulses 
yielded an HIC index for the blackbody of g '^ 4 for 
each interval. In those intervals, TZ was approximately 
constant indicating the emission size of the photosphere 
was constant. Owing to the small number of time bins, 
it is difficult to quantify the evolution of TZ for the single 
pulse GRBs of this study with the coarse time binning 
used, but the fact that the HIC index for the blackbodies 
differs from q ^ 4 and that TZ increases indicates that the 
evolution of the photosphere is very complex. 

5. DISCUSSION 
5.1. The Importance of Fitting with Physical Models 

With the use of physical emissivities in GRB spec- 
troscopy we have shown that the Fermi data are consis- 
tent with synchrotron emission from electrons that have 
not cooled (i.e, slow-cooling spectra) and are inconsistent 
with synchrotron emission from electrons that are cooling 
(fast-cooling). The method enables us to test conclusions 
drawn from the empirical modeling that has dominated 
the field. There is a positive correlation between hard 
alpha and the inability of the data to fit the synchrotron 
model, but the low-energy index of synchrotron seems to 
center ar ound —0.8 rather th an the asymptotic —2/3 as- 
sumed in iPreece et al.l ()1998D . Not only do the fits with 
Band a centering around —0.81 lead to better fits with 
synchrotron, but simulated synchrotron spectra are best 
fit with a Band function having an a sa —0.8. This fact 
alone has very powerful implications for making physical 
interpretations from the value of Band function parame- 
ters. 

Previously, GRB spectra have been successfully fit- 
ted with a thermal -I- non-thermal model , i.e., by us- 
ing a blackbody function combined with a Band func- 
tion. A thermal spectral component has indeed been 
shown to be significant in several cases, foremost in 
GRB 100724B a r id in G RB 110721A (jGuiriec et al.|[20Tll: 
lAxelsson et al. I l2012f ). The Band function in these 
fits is, however, not based on any physical arguments 
but is merely an empirical function that has a pliant 
parametrization. A general problem that arises in this 
type of fitting is that the low-energy slope of the Band 
function (a) and the strength of the blackbody are de- 
generate to varying degrees. When the data is fit by a 
Band function alone, even if an additional component 
really exists in the data at lower energies, it may not 
be identified because the a of the Band function can ac- 
commodate the additional low energy flux by changing 
its slope. 

The slow-cooling synchrotron model we use here is 
more restrictive compared to the Band function. In par- 
ticular, the low energy slope (a) and the curvature of the 
spectrum are predicted by the model. 

The spectrum below the synchrotron vF^, peak is not 
always satisfactorily fit using just the synchrotron model. 
The residuals suggest that an additional component is 
sometimes needed, and we use simulations to show that 
this additional component is statistically significant. Ad- 
ditional components can also be favored in GRBs fit with 
the Band function but we found that the significance of 
the additional component can be greater when added to 
the synchrotron fits than when added to the Band fits. 



We explain this by noting the Band function can accom- 
modate the extra emission using a suitable power-law 
index a but the synchrotron function is more restrictive 
and the additional component is thus required to give 
a reasonable fit. Another essential observation is shown 
in Figure 1101 when using the synchrotron model, the 
temporal evolution of the blackbody flux exhibits well- 
defined pulses and a spectral evolution that is clearly sep- 
arated from the non-thermal emission. This is in contrast 
to the less smooth blackbody flux variations when using 
the Band function as the non-thermal process. This fact 
again reflects that the Band function is less restrictive 
than the synchrotron function and thereby gives rise to 
further scatter in the derived fluxes in the light curves. 
We also argue that these facts suggest that: 

(i) the synchrotron function is a good physical model 
to use 

(ii) the thermal component does exist and 

(iii) multi-component fitting with the Band function can 
be misleading. 

5.2. Alleviating Slow-cooling Synchrotron Problems 

The fact that the non-thermal spectra seem to be 
consistent with slow-cooled synchrotron and not fast- 
cooled synchrotron plac es strong constrai n ts on the emis- 
sion model of GRBs. iGhisellini et al. I ()2000f ) showed 
that trying to reconcile fast-cooling with the obser- 
vations leads to even further problems for the in- 
ternal shock model. In our simple model of fast- 
cooled synchrotron, we neglected the effects of inverse- 
Compton cooling in the Klein-Nishina regime, which 
can significantly alter the value of the low-energy slope 
(jDaigne. Bosnjak. fc Dubus 1 120 lit ). This could make 
some spectra consistent with fast-cooling, but the ef- 
fect on the spectral curvature has not been assessed and 
therefore further investigation is required. Neglecting 
Klein-Nishina effects, the internal shock scenario cannot 
be reconciled with our observations. Additionally, lyyani 
et al. (2013, in prep.) found that for GRB 110721A, the 
standard slow-cooling synchrotron scenario from impul- 
sive energy input such as internal shocks places the non- 
thermal emission region below the photosphere. This 
calculation is possible due to the observation of the black- 
body component, highlighting the importance of multi- 
component fitting. 

Models with co- acceleration via first-orde r and 
second-order Fermi acceleration (jWaxman I 119951 : 
iDermer fc Humill200lD have the ability to balance syn- 
chrotron cooling with stochastic heating keeping 7cooi 
above 7min. The spectrum would resemble a slow-cooled 
synchrotron spectrum. The internal-coll isional- induced 
magn etic reconnection model (ICMART) (jZhang fc Yan I 
120 lit) makes use of magnetic reconnection to accelerate 
electrons to high energies. Electrons are not accelerated 
just once. Instead they are constantly accelerated via 
second-order turbulence. This leads to a balance be- 
tween cooling and acceleration, producing a synchrotron 
spectrum that is in the slow-cooling regime. In addition, 
the ICMART model extends the non-thermal emission 
site far above the photosphere which eliminates the 
problem mentioned above. While simulated spectra from 
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the ICMART model are currently under development 
(Bing Zhang, private communication), none have been 
published to be tested. Still, our analysis suggests that 
such slow-heating scenarios can explain the data much 
better than the current internal shock model. 

5.3. Conclusion 

We have demonstrated that for a set of GRBs we can 
fit a physical, slow-cooling synchrotron model directly to 
the data. Most spectra require a sub-dominant black- 
body with a temperature that places its peak below the 
synchrotron vFi, peak. The temporal evolution of both 
radiative components allows for the assessment of GRB 
jet proprieties without the assumptions required when 
fitting GRB spectra with empirical functions. Several 
parameters in our model can not be constrained by the 
fits, namely 7th, 7min, and B. Detailed assumptions about 
jet properties are required to specify these parameters 
which we will assess in future work. Still, the temporal 
evolution of both synchrotron and blackbody component 
imply that in these GRBs, a photosphere is formed be- 
low a non-thermal emitting region that forms at a larger 
radius than the classical internal shock scenario. The 
electrons in the non-thermal emitting region must un- 
dergo continuous acceleration to produce a slow-cooling 
synchrotron spectrum. These findings help to prune the 
numerous models of GRB emission and provide impetus 
for the further use of physical fitting models. 
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Figure 1. The physical electron spectrum used for fitting synchrotron spectra directly to the data. Here 7th = 10 and 7inin = 30 which 
are the fixed values used for all fits in this study. 




Figure 2. The fast-cooled electron distribution with ')cool = 1 a-^d 7min = 30. A low-energy thermal distribution is excluded in this 
distribution because the synchrotron omission from the power-law form is already too broad for the typical uFv peak of the data and 
including a thermal distribution would smooth it further. 
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Figure 3. The energy flux lightcurvcs of the synchrotron component for the entire sample {black curve). The integration range is from 10 
kcV - 40 MeV for all GRBs except GRB 081224A and GRB 110721A which are from 10 keV - 300 MeV. Superimposed is the slow-cooled 
synchrotron Et (see eq. [5J [red curve) demonstrating the hard to soft evolution of the bursts. 
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Figure 4. The spectral evolution of GRB081224A is an example of the typical evolution observed for the entire sample. The synchrotron 
(from light blue to dark blue) and blackbody (from yellow to red) both evolve from hard to soft peak energies writh time. For this GRB, 
the high-energy power-law corresponding to the electron spectral index does not evolve significantly over the duration of the burst. 



Times: 1.728: 2.752 s 



(U 



3 
O 
O 

CU 
O 



10.0000 
1 .0000 



^ 0.1000 



0.0100 
0.0010 
0.0001 



D 

E 

bo 



2 



■2 




G0_01 + 
NAI_06 O 
NAI_07 
NAI_09 O 
NAI_1 1 '\ 



^1 



^ 



-H 1 — I I I I 



'^^^%m¥f^^^ 










00 1000 

Energy (keV) 



0000 



Figure 5. A time bin of GRB 110721A demonstrating a typical count spectrum of a time bin in the sample. The response has been 
convolved with synchrotron (eqO and a blackbody to produce counts. The residuals from the fit indicate that the model is fitting the 
data well. 
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Figure 6. The distribution of electron indices from the slow-coohng synchrotron fits. Only indices that were constrained are plotted. The 
distribution is broad but centered at S =5 which is much steeper than expected from simple relativistic shock acceleration. 
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Figure 7. The fast-cooled synchrotron fits are poor for nearly all of our sample because none of the spectra have a low-energy index as 
steep as —3/2. Therefore the fast-cooled synchrotron spectrum is too broad around the i/F,^ peak as shown in this example spectrum. 
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Figure 8. An example time bin of GRB 110407A comparing the fitted uF^ spectra of the Band Function [green), slow-cooled synchrotron 
{red), and fast-cooled synchrotron {blue). While the Band function and slow-cooled synchrotron fits resemble each other, the fast-cooled 
synchrotron fit is only able to fit the low-energy part of the spectrum. Because fast-cooled synchrotron has an index of —2/3 below the 
cooling frequency, the fitting engine pushes the value of E» very high to fit the low-energy part of the spectrum resulting in a —3/2 index 
near the uFu peak. The high-energy power-law of the fast-cooling synchrotron spectrum is pushed out of the data energy window. 
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Figure 9. Derived values of the parameter Epeak (obtained using the Band function to fit GRB spectra), versus Et, (obtained using an 
optically-thin non-thermal synchrotron to fit GRB spectra). 
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Figure 10. A subset of flux lightcurvcs illustrating both the temporal structure of the different components and the advantages of 
using a physical model to deconvolve the detector response. The left column contains the lightcurves using synchrotron {blue thick line) 
and blackbody {red thin line) while the right column contains the lightcurves made from using the Band function {green thick line) and 
blackbody {red thin line). The total flux lightcurve {black dotted line) of both approaches are the same. The components have a very 
simple and constrained evolution when using synchrotron as the non-thermal component. This is potentially indicative that synchrotron 
is the actual emission mechanism and the response is being properly deconvolved. In contrast, the lightcurves where the Band function is 
used have large errors and the blackbody does not have a consistent evolution. (See the online version for color) 
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Figure 11. The non-thermal emission of all of the bursts in the sample loosely follow Fe-E, and B»-fluencc relations. For synchrotron, 
the index of this relation should be 1. See Table [5] for the numerical results. 
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Figure 12. The HlC and HFC correlations for the blackbody are separate from those derived from the synchrotron component. This 
adds more evidence for the presence of the component. However, the HIC for the blackbody is not q=4 as expected. 
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Figure 13. The time evolution of kT for four of the GRBs in our sample. GRB110721A is shown w ithout a fit because the coarse time 
binning used did not allow for constraining the fit parameters. However, in lAxelsson et al. I II2012I ). the evolution is shown to follow a 
broken power-law. 
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Figure 14. The evolution of the 7^ parameter is increases with time and shows no relation to the photon flux of the blackbody component. 
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GRB 


Peak Flux (p/s/cm^) 


Duration Analyzed (s) 


Blackbody Component 


LLE data 


GRB 081110A 


20.88 


4.61 


X 


X 


GRB 081224A 


17.11 


18.36 


/ 


/ 


GRB 090719A 


26.52 


30.09 


/ 


X 


GRB 090809B 


18.36 


14.64 


X 


X 


GRB 100707A 


18.77 


22.39 


/ 


X 


GRB 110407A 


15.6 


20.48 


X 


X 


GRB 110721A 


29.82 


12.7 


/ 


/ 


GRB 110920A 


8.08 


238.29 


/ 


X 



Table 1 

The GRBS in our sample. The peak fluxes were taken from the brightest bin of each GRB with a duration determined by Bayesian blocks. 



Time Bin 


Band C-Stat 


Band-I-BB C-Stat 


Synchrotron C-stat 


Synchrotron+BB C-stat 


-0.2-0.2 


453 


450 


485 


455 


0.2-0.8 


374 


362 


695 


364 


0.8-2.4 


427 


405 


2546 


487 


2.4-3.0 


408 


400 


903 


413 


3.0-4.3 


431 


415 


1214 


430 


4.3-5.7 


397 


363 


791 


399 


5.7-7.2 


411 


390 


598 


422 


7.2-12.8 


488 


414 


829 


447 


12.8-22.2 


558 


412 


594 


423 



Table 2 

The time-resolved C-Stat values for GRB100707A show that while the Band function and synchrotron models combined with a blackbody 
function both fit the data well, the non-therr nal fun ctio ns fit t he data very differently when not combined with a blackbody. Specifically, 
where the blackbody is the brightest (Figures |10(c)| and |10(d)| intervals 2, 3, and 4) the Band function alone fits the data acceptably while 

the synchrotron model alone fits the data poorly, mis snows that the flexibility of the Band function can mask the need for the 
blackbody component. The Band-f-blackbody fits actually fit the data better when the blackbody is very bright in this case. This is most 
likely due to the blackbody function (eq. Ilip used is simplified and the actual emission may be broadened due to beaming effects that are 
only important to the fit when the blackbody is bright and the synchrotron fit is used. The Band function makes up for these effects by 
having a harder a. We tested using an exponentially cutoff power-law combined with the synchrotron model and the fits were as good as 
those with the Band function. We will examine the use of a more realistic photosphere model in future work. 



bJO 
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"a; 
=3 



Time Bin 


Band C-Stat 


Band-I-BB C-Stat 


Synchrotron C-stat 


Synchrotron-I-BB C-stat 


Fast C-stat 


Fast-t-BB C-stat 


-0.07-0.08 


640 


640 


673 


673 


709 


709 


0.08-0.48 


690 


690 


704 


704 


1088 


1088 


0.48-1.28 


709 


668 


688 


670 


1654 


957 


1.28-2.78 


887 


761 


838 


770 


1646 


1041 


2.78-3.78 


678 


642 


655 


643 


797 


666 


3.78-5.88 


648 


631 


677 


634 


694 


660 


5.88-7.63 


729 


728 


733 


721 


773 


724 


7.63-12.63 


932 


693 


693 


692 


756 


698 



Table 3 

The C-stat values for ORB 110721A. The significance of the addition of the blackbody is not as large as with GRB 100707A (Table \2i due to the weakness of the blackbody 
component. The fits for fast-cooled synchrotron are included to demonstrate the poor quality fits that are obtained both with fast-cooling synchrotron and fast-cooling synchrotron 

with a blackbody. 
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Time-rcsolvcd Ae 


alysis of Fermi GRBs with Physical Synchrotron Photon 


Models 


Time Bin 


Band a 


Slow-cooled Synchrotron C-stat Fast-cooled Synchrotron C-stat 


■^c-stat 


-5.38-2.82 


-0.9 


523 599 


76 


2.82-3.84 


-0.7 


507 604 


97 


3.84-4.86 


-0.8 


506 596 


90 


4.86-6.91 


-1.0 


534 626 


92 


6.91-9.98 


-1.1 


591 639 


48 


9.98-15.1 


-1.5 


494 494 
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Table 4 

For each time bin of GRB 110407A we examine the C-stat value of synchrotron and fast-cooled synchrotron. This GRB did not have a 
blackbody in its spectrum. While the Band function and slow-cooling synchrotron fit the spectrum well, fast-cooling synchrotron does not 
fit the spectrum unless Band a = —1.5. In this case, the fast-cooled synchrotron peak energy was very unconstrained due to the curvature 

of the data being narrower than the photon model's curvature. 



GRB 


Flux Index 


v2 


$0 


.y2 


GRB 081110A 


2.32±0.4 


0.6 


97±23 


0.4 


GRB 081224A 


1.74±0.1 


1.5 


253±23 


0.3 


GRB 090719A 


1.14±0.07 


0.98 


245±17 


1.2 


GRB 09080B 


1.58±0.05 


8.0 


188±9 


1.0 


GRB 100707A 


1.04±0.02 


1.2 


444±24 


7.3 


GRB 110407A 


1.72±0.20 


0.5 


214±32 


4.2 


GRB 110721A 


1.08±0.03 


14.4 


269±13 


15.4 


GRB 110920A 


1.37±0.06 


0.5 


669±33 


1.2 



Table 5 

Sample correlations for both flux and fluence for the synchrotron component. 



Time Bin 


Band-BB Ac-stat 


Synchrotron-BB i^c-atat 


-0.2-0.2 


3 


30 


0.2-0.8 


12 


331 


0.8-2.4 


22 


2059 


2.4-3.0 


8 


490 


3.0-4.3 


16 


784 


4.3-5.7 


34 


392 


5.7-7.2 


21 


176 


7.2-12.8 


74 


382 


12.8-22.2 


146 


171 



Table 6 

The ^c—stat between the Band function and synchrotron model fits with and without the inclusion of a blackbody for GRB 100707A. 
The blackbody has a significantly larger impact on the fit when included with the synchrotron model. 



Time Bin 


Band-BB l\c-stat 


Synchrotron-BB Ac-stat 


Fast-BB Ac-ataf 


-0.07-0.08 











0.08-0.48 











0.48-1.28 


41 


18 


697 


1.28-2.78 


126 


68 


605 


2.78-3.78 


36 


12 


131 


3.78-5.88 


17 


43 


34 


5.88-7.63 


1 


12 


49 


7.63-12.63 


239 


1 


58 



Table 7 

The i^c-atat values for GRB 110721A tell a different story than GRB 100707A (Table |6j, though both GRBs show a significant 

improvement in the fit when a blackbody is included. Even thought the fast-cooled fits showed extreme improvement with the inclusion of 

a blackbody, the fits arc still poor compared with the slow-cooled model (See Table |3]l. 



GRB 


Flux Index 


v2 


*o 


'^TfA 


GRB 081224A 


2.3 ± 0.3 


1.4 


121 ± 13 


9 


GRB 090719A 


2.8 ± 0.4 


2.3 


232 ± 23 


3 


GRB 100707A 


2.2 ± 0.1 


17.4 


319 ± 8 


28 


GRB 110721A 


1.3 ± 0.2 


1.8 


43 ± 3 


5 


GRB 110920A 


2.0 ± 0.1 


0.9 


1147 ± 21.7 


4 



Table 8 

For the subset of bursts that have a strong blackbody component wc compute the flux and fluence correlation for the blackbody. 
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GRB 


Fbb/Fsyn 


1^* Decay Index 


2"'* Decay Index 


v2 


GRB 081224A 


0.3 


-0.6 ± 0.07 


-20 ± 243 


0.5 


GRB 090719A 


0.4 


-0.1 ± 0.05 


-2.0 ± 0.7 


11.3 


GRB 100707A 


0.5 


-0.4 ± 822749 


-0.8 ± 0.03 


22.7 


GRB 110920A 


0.8 


-0.3 ± 0.03 


-0.9 ± 0.04 


2.0 



Table 9 

The evolution of the blackbody follows a broken power-law. However, the coarse time bins recovered by the Bayesian blocks algorithm 

make it difficult to constrain the decay indices. 



